Entropic forces generated by grafted semifiexible polymers 
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The entropic force exerted by the Brownian fluctuations of a grafted semifiexible polymer upon a 
rigid smooth wall are calculated both analytically and by Monte Carlo simulations. Such forces are 
thought to play an important role for several cellular phenomena, in particular, the physics of actin- 
polymerization-driven cell motility and movement of bacteria like Listeria. In the stiff limit, where 
the persistence length of the polymer is larger than its contour length, we find that the entropic 
force shows scaling behavior. We identify the characteristic length scales and the explicit form of 
the scaling functions. In certain asymptotic regimes we give simple analytical expressions which 
describe the full results to a very high numerical accuracy. Depending on the constraints imposed 
on the transverse fluctuations of the filament there are characteristic differences in the functional 
form of the entropic forces. In a two-dimensional geometry, the entropic force exhibits a marked 
peak. 

PACS numbers: 05.20.-y, 36.20.-r, 87.15.-v 



I. INTRODUCTION 

In a cellular environment soft objects like membranes 
and polymers are subject to Brownian motion. As a re- 
sult there are interactions between them which are en- 
tropic in origin, i.e. a consequence of constraints imposed 
on the Brownian fluctuations. For example, two parallel 
membranes repel each other entropically with a poten- 
tial that falls off like a power law in the distance between 
them [IJ. Similarly, thermally fluctuating biopolymers 
like F-Actin and microtubules may exert entropic forces 
on membranes or some other obstacles; for an illustra- 
tion see FigQ] Though due to the same thermal fluctu- 




FIG. 1: A smooth hard wall constrains the Brownian fluctua- 
tions of a grafted semifiexible polymer resulting in an entropic 
force on the wall. 



ations such forces have to be distinguished from forces 
obtained by pulling on a biopolymer 0- 01 It will turn 
out that the force-distance curves of these two cases have 
no resemblance at all in a regime where thermal fluctu- 
ations play a role, which is generically the case for all 
cytoskeletal filaments. Both types of forces are thought 



to play a prominent role in cell motility and movement 
of pathogens like listeria monocytogenes, that propel it- 
self through the cytoplasm of infected cells by construct- 
ing behind it a polymerized tail of cross-linked actin fil- 
aments 4]. Similarly, in a crawling cell, the force gen- 
erated from the growth of a collection of actin fibers is 
responsible for the protrusion of cell membrane, which 
are known as lamellipodia, filopodia, or microspikes ac- 
cording to their shapes [5j- It seems that quite gener- 
ally polymerizing networks of actin filaments are capable 
of exerting significant mechanical force, which are used 
by eukaryotic cells and their prokaryotic pathogens to 
change shape or move. One type of force is generated by 
fluctuating filaments at the leading edge of the network. 
The length of the thermally fluctuating parts of these 
polymers are typically 200 ~ 300nm, which is very short 
compared to their persistence length £ p sa 17/im 6], such 
that an analysis which considers these filaments as stiff 
seems appropriate. 

In this paper we will not enter into the debate on the 
particular force generating mechanism responsible for all 
these different types of cell motility, but rather give a 
detailed analysis of the entropic forces which fluctuat- 
ing stiff polymers exert on rigid walls. This may serve 
as important input for future molecular models of force 
generation in cellular systems. 

Consider a semifiexible polymer with contour length L 
and persistence length £ p with one end fixed both in po- 
sition and orientation to some rigid support (see Fig^l, 
e.g. the dense part of an actin gel. We choose coor- 
dinates such that the grafted end is at the origin with 
the tangent fixed parallel to the z-axis. Consider a rigid, 
smooth wall orthogonal to the a>z-plane at a distance £ 
from the origin. Let i? be the angle between the z-axis 
and the normal h of the wall. If C is small enough, the 
wall will constrain the Brownian fluctuations of the poly- 
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mer leading to an increase in free energy with respect to 
the unconstrained polymer. On time scales larger than 
the equilibration time of the grafted polymer this results 
in an average force / exerted on the wall. Our goal is to 
calculate how the entropic force / depends on £ and ■& 
and the contour length L and the persistence length £ p 
of the polymer. 

We will proceed as follows. The following Section [H] 
serves to introduce and discuss the various types of ther- 
modynamic forces which can be generated by fluctuating 
scmiflexible polymers. We will arrive at the conclusion 
that the entropic forces discussed above are closely re- 
lated to the probability distribution of the free end of 
the clamped polymer. In Section ITTT1 we start our analy- 
sis of entropic forces with a polymer grafted perpendicu- 
lar to the wall. This chapter contains a definition of the 
wormlike chain model and the basic idea of our analytical 
calculations, which starting from the tip distribution cal- 
culates the restricted free energy and the entropic force. 
The analysis is complemented by Monte Carlo (MC) sim- 
ulations, which both show the range of validity of the 
analytical results and the crossover from semiflexible to 
Gaussian chains. Details of the calculations are deferred 
to the Appendices EI IEI an d 10 Section IIVI treats the 
technically more complicated case of a polymer inclined 
at an angle d with respect to the wall. Here we obtain the 
entropic forces analytically up to the numerical evalua- 
tion of some integrals. For some asymptotic cases explicit 
analytical formula are again obtained. The MC simula- 
tions in this chapter are restricted to a parameter range 
which is close to the stiff limit, and mainly serve the pur- 
pose to define the range of applicability of the anayltical 
results. Finally, in the conclusion we give a discussion of 
our main results. 



II. ENTROPIC FORCES AND PROBABILITY 
DENSITIES 

According to the wormlike chain model 0, 0] , the elas- 
tic energy of a given configuration r(s), parameterized in 
terms of the arc length s £ [0, L], is given by 

Here t(s) = dr(s)/ds = r(s) is the local tangent to the 
contour r(s) and £ p = n/k^T is the persistence length 
with k the polymer's bending modulus, and f3 = 1//cbT. 
As the polymer is considered to be inextensible, we have 
|t(s)| = 1 for all s, i.e. the tangent vectors are restricted 
to the unit sphere. 

In a cellular environment biopolymers are flexed by 
Brownian motion, i.e. they exhibit thermal fluctuations 
in their shape. This mere fact makes for a rich me- 
chanic response genuinly different from its classical ana- 
logue, a rigid beam. Consider a polymer whose posi- 
tion (not its orientation) is fixed at one end and one 



is pulling on its other end, a typical situation encoun- 
tered in an experiment using optical or magnetic tweez- 
ers. Then there is no unique force-distance relation. It 
actually matters whether one pulls at constant force / 
and measures the resulting average distance (r)(f) or 
vice versa. Results for the constant force ensemble are 
thoroughly discussed in Ref.Q- In a constant distance 
ensemble the probability density distribution of the end- 
to-end distance P(r) provides the necessary information 
0. It defines a free energy F(r) = —k^T In P(r) from 
which the average force may be derived by differentiation, 
(/)(r) = -dF(r)/dr 

Here we are interested in the force a fluctuating fila- 
ment exerts on a rigid obstacle which is fixed in its posi- 
tion • The polymer's end facing the obstacle is consid- 
ered as free to fluctuate and only its proximal end is fixed 
in position and orientation; see Fig^ Since there are no 
direct forces between polymer and obstacle the force ex- 
erted on the wall is solely due to the steric constraints 
imposed on the filament. This suggests to use the term 
"entropic forces" , frequently used in analogous physical 
situations [T^. However, this should not leave the reader 
with the wrong impression that there are different phys- 
ical origins for entropic forces and those discussed in the 
preceeding paragraph. It is merely the type of "boundary 
condition" imposed on the thermal fluctuations which 
leads to their (drastically) different character. 

For getting acquainted with the problem let us consider 
the simplest case, a grafted polymer whose one end and 
tangent is fixed such that it is oriented perpendicular to 
a smooth wall (Fig^with d — 0). The presence of the 
wall allows only for those polymer configurations which 
are entirely in the halfspace to the left of the wall. Since 
we are mostly interested in stiff polymers (which have 
a low probability for back-turns) this restriction may be 
approximated as a constraint solely on the position of the 
polymer tip facing the wall, r z (L) < £; later in Section 
IIV CI we will show some simulation data going beyond 
this approximation. 

To derive the average force acting on the wall we con- 
sider a wall potential [/(£ — r z (L)) for the free polymer 
tip, which at the end of the calculation will be reduced to 
a hard wall potential. For now picture a steep potential 
which rises rapidly for r z (Li) — > Q. Then, the ensemble 
average for the force the polymer tip excerts perpendic- 
ular to the wall reads 

Here the partition sum 

Zy(C) = Jv[r( s )}e-^ H+U ^ (3) 

is a path integral over all polymer configurations compat- 
ible with the boundary conditions imposed on the distal 
and free end of the grafted polymer, where the measure is 
taken such that the partition sum without a constraining 
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wall (U = 0) is normalized to 1. This is now a thermo- 
dynamic force. In an actual experiments it is obtained 
by a time average with an averaging time much larger 
than the equilibration time for the grafted polymer. This 
force would also be measured in an experiment where a 
large number of independent and identical polymers push 
against the same wall. 

Since the wall potential depends only on the difference 
between the position of the polymer tip and the wall we 
may rewrite the entropic force in Eq. (J2J) as 



(/||)(C) = fc B T^lnZ||(C). 



(4) 



Upon definig a free energy of the confined polymer as 



^11 (C) = -fc B TmZ||(C) 



(5) 



the entropic forces again reads as a spatial derivative of 
a free energy 



(/||XC) = -^ii(C) 



(6) 



The physical interpretation of this free energy becomes 
clear as one goes to the hard wall limit. Then, the parti- 
tion function reduces to 



Z|,(0 = J V[r( S )]Q(C-r z (L))e 
= : (0(C-r,(i)))o, 



(7) 



where the subscript indicates that the average is now 
taken with respect to the bending Hamiltonian only. The 
O-function, defined such that 8(x) = 1 for x > and 
zero elsewhere, indicates that only those configurations 
are counted with the position of the polymer tip to the 
left of the wall. Hence, as for the fixed distance ensemble 
in a pulling experiment, the free energy results from a 
quantity measuring the number of configurations obey- 
ing the imposed constraint, where each configuration is 
weighted by a Boltzmann factor for the bending energy. 
It is useful to rewrite the partition function as 

Z„(0 - f dze>(t-z)(S(z-r z (L))) 
K 

dzP\\(z), (8) 

-L 

where P\\(z) = (8(z — r z (L)))o is the probability density 
to find the z-coordinatc of the polymer's free end at z, 
irrespective of its transverse coordinates. It identifies the 
restricted partition sum as the cummulative distribution 
function corresponding to the probability density P\\{z). 
One may then write the entropic force in the alternative 
form 



</||}(C) = k B T 



P\\(C) 



(9) 



Upon multiplying this formula by d£ it may be inter- 
preted as follows. The work done on the wall upon 
displacing it by an infinitesimal distance dC equals the 
thermal energy scale k B T times a conditional probability 
-FWt(0^C = P\\ (C) t ^C/-2'|| (C); which measures the prob- 
ability that the position of the polymer tip is within a 
distance d£ from the wall given that the polymer is in 
the left halfspace. 

Since the probability density for the position of the 
polymer tip -P(x, z) is actually a function of the position 
perpendicular and transverse to the wall, Eq. J^J immedi- 
ately suggests that one could define a local entropic pres- 
sure. Indeed, upon generalizing the above arguments one 
may write 



P(x,0 



*q Jv[r]^5{x-r ± {L))e-n*-W 
(0(C-r,(L))a(x-rx(L))>o 



3i(0 

k B T d 

MO 



knT 



(10) 



for the entropic pressure, i.e. the force per unit area 
exerted locally at x on the wall. Again, the entropic force 
is given by the thermal energy scale times a conditional 
probability density, which now measures the probability 
of finding the polymer tip at a particular site x on the 
wall conditioned on the polymer configuration being to 
the left of the wall. Pictorially, one may say that the local 
pressure is given by k B T times the number of "collisions" 
of the polymer with the wall per unit area, a reasoning 
which is frequently used in scaling analyses. 

The total force is, of course, obtained by integrating 
over this local pressure, (/ii)(C) = J dxp(x, Q. In addi- 
tion, one may now also define an entropic torque as has 
recently been done for a rigid rod facing a planar wall 
|l3| : we leave this issue for future investigations. 

Generalizing the above ideas suggests to introduce an 
effective local free energy per unit area as 



JP(x,C)--fc B T 



(ii) 



which is useful in applications where the obstacle is ac- 
tually not rigid but soft with some internal elasticity, 
e.g. a membrane, whose dynamics is much slower than 
the equilibration time of the polymer. Then the elas- 
tic energy describing membrane bending and the above 
effective free energy may just be added to describe the 
combined system. Of course, such a description fails if 
time scales for the dynamics of both soft objects are com- 
parable. 

Our main conclusion in this section is that entropic 
forces generated by a grafted stiff polymer can be re- 
duced to the calculation of the probability distribution 
of the polymer tip. For a polymer constrained to two 
dimensions this distribution function has been found to 
show quite interesting behavior such as bimodality in the 
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transverse displacement of the free end 0]. This pro- 
nounced feature of the distribution function has recently 
been rationalized upon exploiting an interesting analogy 
to a random walker in shear flow |l5| . 



III. POLYMER ORTHOGONAL TO A WALL 

In this section we are going to calculate the entropic 
force generated by a grafted polymer whose orientation 
is on average perpendicular to the wall. It illustrates the 
basic idea of our analytical calculations for the simplest 
geometry. 



A. Weakly bending limit: mode analysis 

In evaluating the distribution function analytically we 
restrict ourselves to the limit of a weakly bending filament. 
In other words, we consider the persistence length £ p to 
be large enough compared to the total contour length 
L, such that the statistical weight of configurations with 
small sharp bends will be negligible. The key small di- 
mensionless quantity will be the stiffness parameter 



e = L/£ p 



(12) 



and we will refer to the weakly bending limit also as the 
stiff limit. 

For small e, the transverse components, t x (s) and 
t y (s), of the tangent vector t(s) will be small for all 
s. While the condition |t(s)| = 1 would suggest a pa- 
rameterization of t(s) in terms of polar coordinates or 
Euler angles, for reasons that will become apparent later 
we want to have independent variables that appear in 
a symmetric way in the integrand. Thus we choose to 
parameterize t by 



(13) 



1 




where we dropped all arguments s for brevity; the gen- 
eralization to d spatial dimensions is obvious. 

The boundary conditions at the ends of the polymer 
are 

t(0) = (0,0, 1) T (clamped end) , (14a) 
t(L) = (0,0, 0) T (free end). (14b) 

This translates into a(0) = (a x (0), a y (0)) T = (0, 0) T and 
a(L) = (d x (L),dy(L)) T — (0,0) T . We thus can choose a 
Fourier representation, or in other words a normal mode 
decomposition 



oo 



(15) 



with eigenvalues 



(2k- 1) 



(16) 



and Fourier (normal mode) amplitudes 



a x .k 



L 



ds 



a x (s) sin (^kjj , (17) 



and similar for a y (s). To second order in the Fourier 
amplitudes the location of the end-point along the 2-axis 
reads 



r z (L) = / dst z (s) 
Jo 







ds [a 2 x (s) + a 2 y (s)] 



L--Y 
4 ^ 

k=l 



a x,k + a y,k\ 



(18) 



Similarly, we find for the Hamiltonian to second order 
/^«|ff>I[< fe + < fc ]. (19) 



k=l 



B. Moment generating function 

For calculating the probability density function P|| (z) 
we follow a procedure outlined in Ref . [|J and consider the 
moment generating function 



/ e -/(L-r.(L)) )() 



-L 
2L 



dze- i{L - z) P\\{z) 



/ dpe-^P^L-p). 
Jo 



(20) 



Note that thermal averages have to be evaluated using 
the bare elastic free energy, Eq. Since for stiff chains 
configurations with large values for the stored length 
("compression") p = L — z are rather unlikely, we can 
extend the upper boundary of the integral in the last 
line of the preceding equation to infinity. This allows us 
to write the moment generating function as the Laplace 
transform of the distribution function Pm (z) 



dpe- fp P«(L- p). 



(21) 



For / = the latter equation reduces to the normaliza- 
tion condition of the probability density function P\\ (z) 
such that P||(0) = 1. 

Combining Ea s . 1 1 81 [HH and I2TH the moment generating 
function can be put into the following path integral form 



k=l 



J V[ a (s)} exp j-i £ds [£ p a 2 + fa 2 ] j (22) 
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with the boundary conditions given by Eq. I|14|), This 
path integral is easily evaluated upon using the Fourier 
representation of the transverse tangent fields, Eq. (|15Jl . 
and noting that to harmonic order fluctuations in all 
transverse directions are statistically independent. We 
find in d spatial dimensions 



k=l 



naa-k 
^exp 



A k l P 

L 



fL 




n 



2 n, -(d-l)/2 



(23) 



where the normalization factor J\f of the path integral 
was chosen such that Vn (0) = 1. If / € M + the product 
may be rewritten as |l7j 




(24) 



Note that the moment generating function, which also 
depends on the length scales L and £ p , has the scaling 
form 



(25) 



where we have defined the characteristic longitudinal 
length scale 



L 



(26) 



The formulas Eq. and Eq. are the basis for all 
subsequent calculations in this section, which are basi- 
cally different forms of performing the inverse Laplace 
transform. 

For future reference and comparison with the entropic 
forces we close this subsection with a discussion of the 
force-extension relation in the fixed force ensemble. It 
simply follows as the first moment of the moment gener- 
ating function 



= Hi — 



IT 

df 



L(d-1) tanhy/L^ 



4£ n 



(27) 



where / is the external force in units of the thermal en- 
ergy k^T . In the limit of small external forces this re- 
duces to 



:(L))f = L 



1 - 



d- 1 L d-1 



12 



(28) 



which identifies Lh as 4/(d — 1) times the equilib- 
rium stored length due to thermal fluctuations. We 
also recover the effective linear spring coefficient fc|| = 



12k 2 J (d — l)kBTL 4 , which was previously calculated in 
Ref. [18j . For strong stretching forces the extension satu- 
rates asymptotically as 



(r z (L)) f =L 



L (d-1) 



4 W7^i 



(29) 



In the limit of large compressional forces the weakly 
bending rod approximation breaks down and one has to 
use different approaches to evaluate the force-extension 
relation [Io|. 



C. Probability density for the position of the 
polymer tip: analytical and MC results in 3d 

We now return to the distribution function and the 
resulting entropic forces. Upon performing the inverse 
Laplace transform one gets (for details of the calculations 
see Appendix I A 1|) 



L 



£(-l) fe+1 A fc exp 



L 



L 



(30) 



Inspection of Eq. I|30|l immediately tells us that it can 
be written in scaling form 



P«(z,L,£ p )=L^Pn(p), 



(31) 



where we have made the dependence of the probability 
density on L and £ p explicit and introduced the scaling 
variable 



L 



P 



(32) 



measuring the compression of the filament in units of L|| . 
This implies that data for the probability density of the 
polymer tip can be rescaled to fall on a scaling function 
P|| (p), shown as the solid curve in Fig|3 Of course, since 
the analytical calculations are based on the mode analysis 
in the weakly bending limit such a universal scaling curve 
is obtained only for small enough stiffness parameters e. 

The probability density is strongly peaked towards full 
stretching, p — > 0, and falls off exponentially for large p, 
such that for p > 0.3 



1 2~ 
7T exp [ — — 7T p 



(33) 



is already an excellent approximation. The scries expan- 
sion given by Eq. I|30|l converges well for all values of 
z well below L, but its convergence properties become 
increasingly worse if z approaches L. As detailed in Ap- 
pendix one may also derive an alternative series rep- 
resentation of the tip distribution function which con- 
verges well close to full stretching 



21 + 1 



' 7T/5 3 



exp 



1)21 



(34) 
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0.2 0.4 0.6 0.8 1 

p=(L-z)/L|| 

FIG. 2: Scaling function Pit (p) (solid line) in 3d for the prob- 
ability density to find the free end of a grafted semiflexible 
polymer in a plane defined by r z (L) — z or equivalently with 
a reduced stored length p. For comparison MC data are given 
for a series of stiffness parameters e = L/£ p indicated in the 
graph. Deviations from the scaling curve in the stiff limit 
become significant for e = 0.5 and larger. 



Already the first term of Eq. 



'irp 3 



exp 



1 

~4~p 



(35) 



gives an excellent fit for p < 0.3. In particular, it captures 
the main feature of the distribution function, namely its 
maximum close to full stretching. The same approximate 
expression may also be obtained by evaluating the inverse 
Laplace transform using the method of steepest descent; 
see Appendix^] The asymptotic results given in Eq. 
and Eq. (|33() taken together give a representation of the 
scaling curve to a very high numerical accuracy. They are 
the analogues of the results found in Ref.@ for a freely 
fluctuating filament; see also Ref.[To|. 

The MC data shown in Fig|21have been obtained by us- 
ing a standard algorithm for a discretized wormlike chain, 
similar to the one described in Ref.0- As expected, the 
MC results agree very well with the analytical calcula- 
tions for small values of e. From Fig we can read off 
that the asymptotic stiff scaling regime remains valid up 
to stiffness parameters e « 0.1; even for e = 0.5 the 
shape of the scaling function resembles the MC data quite 
closely. As the polymer becomes more flexible the shape 
asymptotically becomes Gaussian; for e = 3 a skew is 
still noticeable. Note that in the parameter range given 
in FigI2 the width of the rescaled probability densities 
stays approximately the same and is hence well charac- 
terized by the longitudinal scale L». 



D. Confinement free energy and entropic forces: 3d 

Now we are in a position to calculate the restricted 
partition sum (cummulative probability distribution) 

Zu (£) = J_ L dzP\\ (z) by (formally) integrating the series 
expansion Eq. (|30(l term by term. This gives 



■ZiiK) 



dz P|i (z) 



OO 

l-2^(-l) fe + 1 A,T 1 (l 



k=l 



2E(- 1 ) fc+1 V e 

fe=l 



-\l{L-0/L\\ 



(36) 



where in the first line we used the normalization of P|| (z) 
and in the final line the identity |l6j | 



k=l 



■1) 



fc+1 



1 



2k- 1 



7T 

4 



(37) 



The series expansion in Eq. (|36|l converges well for all 
values of £ well below L. Alternatively, one may start 
from Eq. I|34|) and derive 



Z|,(0 = l + 2£(-l) fe erfc 



Afe/TT 



k=i 



y/(L-Q/L\\ 



,(38) 



which is well behaved for ( close to L, and dominated by 
its first term. A second method to obtain Eq. (|38|l can 
be found in Appendix [0 

/From both series expansions it is evident that the re- 
stricted partition sum has the scaling property 



Z||(C,L,^ P ) = Z||(r7), 
where we have introduced the scaling variable 

- L-C 



V 



(39) 



(40) 



which measures the minimal stored length (compression) 
t] = L — ( of the filament in units of L». The confine- 
ment free energy, F\\{fi) = —k^T hi Z\\(ij) , corresponding 
to this partition function is shown in FigEJ Again, the 
universal scaling function describes the MC data well for 
e < 0.1. Note that for all values of fj and the stiffness 
parameter e the free energy is convex. This will turn out 
to be an important feature which distinguishes the 3d 
and 2d case. 

Upon using Eq. ^ for the entropic force we find 



/ii (0 



k B TPi\(y) 



(41) 



which immediately shows its scaling behavior and iden- 
tifies ksT/L^ as the characteristic force scale. It is up to 
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0.15 0.2 0.25 

rj = (L-Q/L« 



0.4 



FIG. 3: Confinement free energy JFn (fj) of a grafted polymer 
constrained by a rigid wall in 3d as a function of the reduced 
minimal stored length fj = (L — Q/Lu. The solid line gives 
the scaling function obtained in the limit of a weakly bending 
rod. Symbols represent MC data for different values of the 
stiffness parameter e as indicated in the graph. 



a prefactor identical to the critical force 
tt 2 k tt 2 k n T 



fc 



4L 2 



4 L 



(42) 



for the buckling instability of a classical Euler-Bernoulli 
beam 19] . It suggest to rewrite the entropic force as 



/ll(C,£,*p) = /c/n(»7), 

with the scaling function 



f\\iv) 



4 P\\(fj) 

ir 2 



(43) 



(44) 



The analytical result for the scaling function f\\(fj), 
shown as the solid curve in FigQJ has several charac- 
teristic features. First of all, it is always monotonically 
increasing since the free energy is convex. For fj ^ 0.4, 
the scaling function is /|| » 1 corresponding to / c , 
i.e. a vanishing contribution of thermal fluctuations to 
the force. For smaller fj, corresponding to larger distances 
Q between the wall and the grafted end of the polymer, 
fluctuations reduce the force exerted on the wall by effec- 
tively shortening the polymer. For C, — > L (resp. fj = 0), 
the probability of the polymer to contact the wall be- 
comes smaller and smaller until finally for £ = L only 
one configuration, namely the completely straight one, 
has r z (L) = L. Hence the force must vanish for all ( > L 
(resp. fj < 0). 

We have learned already in Section llll Cl that there are 
excellent approximations to the scaling function for the 
probability density of the free polymer end for small val- 
ues of the reduced stored length, Eq. l|35|l . In the same 




FIG. 4: Scaling function f« (fj) for the entropic force exerted 
on a wall at a distance £ from the grafted end as a function 
of the scaling variable fj = (L — Q/Ln. MC data for different 
stiffness parameters e as indicated in the graph. 



way, the first term of Eq. 1|38|) is an excellent approxima- 
tion to the infinite series for fj ^ 0.2. Thus we may write 
for the scaling function of the entropic force 



/|f(*?) 



4e -l/4"? 



7r 5/2^3/2( 1 _2erfc(l/2 V ^)) 



(45) 



which already describes most of the nontrivial shape of 
the scaling function. For fj ^ 0.2 it suffices to high accu- 
racy to use the first two terms of Eq. I|36|) , which gives 



1 - 3e~ 27r 



1 



Ip-27T 2 77 



(46) 



Upon inspection of Eq. I)41|l one may interpret the func- 
tional form of the entropic force as due to two effects. 
In the numerator we have the probability density for the 
position of the free end at the wall. This function shows 
a pronounced peak as one decreases the distance C (resp. 
increases the scaling variable fj). At the same time, the 
denominator, the cummulative distribution function, de- 
creases with decreasing £. It is now a matter of how fast 
these changes occur what the ensuing shape of the scal- 
ing function for the entropic force will be. In the present 
case of a polymer in 3d the decrease in the cummulative 
distribution function seems to be fast enough to com- 
pensate the maxiumum in the probability density of the 
free polymer end such that the entropic force becomes a 
monotonically increasing function of fj. 

^From Fig0] one observes that the universal scaling 
curve is a lower bound to the MC data for all values of 
the stiffness parameter e. For fixed e the entropic force 
always increases monotonically with increasing compres- 
sion; for intermediate values e ~ 2.5 there is a pro- 
nounced change in curvature at fj sa 0.25. For strong 
compression the results asymptote to the mechanical 
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limit (fceT = 0). This limit is not correctly reproduced 
within the harmonic approximation which gives 



/mech(C) = fcQ(L - C) , 



(47) 



whereas the exact force-extension curve is a monotonous 
function in £ that is somewhat larger than f c for £ < L 
and tends to f c for £ — > L. 

One might finally ask, whether these entropic forces 
/|| (£) are related to the force extension relation discussed 
in section lrTTBl (r z (L)) f - (r z (L)) = k^\f + 0{f 2 ) with 

feii = 6k 2 /kftTL 4 [l8j. Rewriting these linear response 
result in scaling form we find, 



fc 



24 



(48) 



Comparing this with FigQJ we see that the linear re- 
sponse result does not contain any information about the 
situation under investigation here. To the contrary, the 
initial rise of the force when £ becomes slightly smaller 
than L is highly nonlinear (see Eq. (@5jl). 




FIG. 5: Probability density Pit (p) of the free end of a grafted 
semiflexible polymer in 2d as a function of p (solid curve). 
MC data for different stiffness parameters as indicated in the 
graph. 



E. Distribution function and entropic forces: 2d 

Analogous to the previous section the tip distribution 
function of a polymer confined to 2d, e.g. by two parallel 
glass plates, obeys a scaling law in the stiff limit 

P\(z,L,e v ) = L^ 1 P\(p). (49) 

The scaling function may again be represented in terms 
of series expansions (see Appendix IA 2|) . A series which 
converges well for small values of p reads 



_i 

2 

I 



21 + 



2ttp 3 / 2 



exp 



(50) 



for an explicit formula for the binomial coefficient in 
Eq. H5U|) see Eq. I|A17|I . The scaling function, shown as 
the solid curve in Fig|51 has an overall shape which is 
quite similar to 3d with a pronounced maximum close 
to full stretching. The series approximations may again 
give useful approximate expressions for the shape. In the 
proximity of full stretching, the series given by Eq. I|5(J|) 
converges very fast such that already the first term 



if (p) = 



1 



■ exp 



1 



(51) 



^8np 3 / 2 

is an excellent approximation for the whole series at least 
for p < 0.3. As in 3d, a saddle point approximation 
also gives Eq. 1)51(1 (see Appendix [El. Alternatively, as 
shown in Appendix IA 21 one may derive a series expan- 
sion which converges well in the strong compression limit; 
see Eq. (|A11|) . For p g2 0.3 it suffices to use the first term 
of this sum only which reads 



,-7T 2 p/4 



2^2 

2.5 e" 21 ^/ 16 



l + 1.5e- 5 ^ 16 + 2e- 12 ^/ 16 



3e 



-32tt 2 p/16 



(52) 



Upon increasing the stiffness parameter the rescaled 
probability distribution deviates from the scaling func- 
tion in the semiflexible limit and approaches a Gaus- 
sian. In contrast to 3d there is an intermediate parameter 
regime in the stiffness parameter where P11 (p) exhibits a 
marked shoulder. This feature of the distribution func- 
tion has recently been identified and explained in terms 
of an interesting analogy with the physics of a random 
walker in shear flow [15j. 

Upon integrating Eq. I(50|l from — L to £ one obtains 
for the restricted partition sum 

? «\ 1 AT^ (-l) fc (2fc-l)!L ^A 2fc+1 ^ 
2||(0=1-V2L ^ er H^J ' ( } 



fc=0 



with the same scaling variable 77 as in the previous sec- 
tion. Similarly, using Eq. (|Alip gives 

-.00 s 

z \\ (0 - E(-!) fe E ^/4 e " AL+ * /4 " ■ (54) 

fc=0 2=4 

Hence, as in 3d, one finds for the free energy 

^U(C-Mp) = -fe B Tlni||(j7) (55) 
and the entropic force 

/||(C,L,£ p )=/ c /||(^) (56) 
with the scaling function 

4 Z[\{fl) 



tt 2 Z {[ (rj) 



(57) 



and f c = 7t 2 k/4L 2 ; see the solid curves in Fig|S|and Fig[7| 
for a plot of the scaling functions for the free energy and 
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?j = (L-C)/i|| 

FIG. 6: Free energy of a grafted polymer whose tip is con- 
strained by a rigid wall in 2d. The solid line gives the universal 
scaling function in the stiff limit. MC data are given by the 
symbols for different values of the stiffness parameter e as 
indicated in the graph. 



obtains for the entropic force 

- _ y/2e- 1/16i > 
^ ™ ~ ^5/2^3/2 p _ 72erfc[l/4v^]) ' 

This has the same functional form as the corresponding 
expression in 3d, Eq. (|45|) . but differs in some numerical 
factors. These differences can all be traced back to the 
strength a of the essential singularity of the tip distribu- 
tion function close to full stretching, P\\{p) oc exp(— a/p); 
compare Eq. (|35|) with Eq. (|51|1 . One may interpret this 
strength as a kind of phase space factor counting how 
fast the number of polymer configurations decreases as 
one approaches full stretching. It clearly shows that the 
maximum of the entropic force in 2d is of purely geomet- 
ric origin. As an interesting consequence of this maxi- 
mum one should note, that for most values of the reduced 
stored length i) the entropic force exceeds the purely me- 
chanical force given by the Euler buckling force. 

IV. GRAFTED POLYMER AT AN OBLIQUE 
ANGLE TO THE WALL 




V=(L-0/h 

FIG. 7: Scaling function for the entropic force which a grafted 
polymer exerts on a rigid wall in 2d as a function of the re- 
duced stored length fj — {L — (,) / L\\. The solid line gives the 
universal scaling function in the stiff limit. MC data are given 
by the symbols for different values of the stiffness parameter 
e as indicated in the graph. 



entropic force, respectively. The key difference between 
the results in 2d and 3d is that the effective free energy 
exhibits a change in curvature at i) w 0.05 and as a result 
a pronounced peak in the entropic force. The peak is 
a pretty robust feature of the distribution function and 
vanishes only for very large values of e ~ 5. 

In order to understand the physical origin of this peak 
it suffices to consider small values of r). Then, using only 
the leading term of the series expansion Eq. I|53|) . one 



The generic situation one encounters in a cellular sys- 
tem is that the polymer is inclined with respect to a 
membrane. Then we have to ask how the force derived 
above changes when the graft of the polymer is not or- 
thogonal to the constraining wall but at some oblique 
angle it/ 2 — i9; see FigEl Since the presence of the wall 
restricts the position of the polymer tip to 

r z (L) cos i3 + r x (L)smi3 < ( (59) 

one has to evaluate the restricted partition sum 

Z(C,tf) = (e[C-r z (L)costf-r a (L)sintf]) 

= JdxdzPfrzmC-zcctt-xdaq (60) 

to find the entropic force. 

A. Probability distribution function of the tip 

This calculation requires the knowledge of the joint 
probability density of the tip 

P(x, z) := (5[r x (L) - x] 5[r z {L) - z]) . (61) 

In Section IIIII we have already analyzed the reduced 
distribution function Pm (z) and found that its width is 
characterized by the scale Lh = L 2 /£ p . Similarly, one 
can find an explicit expression for P±(x) in harmonic 
approximation, where 

r x (L) w y2a x ,k / ds s'm(X k s/L) 
k=i Jo 

OO 

= L^AfcVfc, (62) 
fe=l 
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FIG. 8: A smooth hard wall at some oblique angle 7r/2 — $ 
constrains the configurations of a stiff polymer grafted parallel 
to the z-axis. 



tribution 



Px(x) 



2ttL 



-e 2 



(64) 



where we have defined the characteristic transverse length 
scale 



L x = JL 3 /U 



(65) 



Together with Lii these are the two length scales are char- 
acterizing the width of the joint distribution function. 
This suggests to write the joint distribution function as 



P(x, z. L, £p) — 



1 



LnL 



-P(x,p) 



(66) 



and thus 



Px(x) = f ^eWexp 



dq 
2^ 



e iqx exp 



i p fc=i 



(63) 



With J2kLi (2fc-i)4 = M E3> *h* s S ives a Gaussian dis- 



in terms of dimensionless variables 

x = x/L± , 

P = ( L - z )/ L \\ 



(67) 
(68) 



An explicit form of the joint distribution function can 
be calculated to harmonic order. For simplicity, we start 
with a polymer fluctuating only in the x-z-plane [d = 2). 
Then 



dq z dq x 



I Hit 

dq z dq x 
2tt 2tt 

dq z dq x 
2tt 2tt 



^1^3 (L — z)—iq x x 



11 (exp 



4 



2 

-1 0C 



jqx{L—z)—iq x x 



n 



■ exp 



3^1 



\i{\l + iq z L\\) 



1 


L 


Li|| 




1 


I 


ill 



nvi+i5 



-2 ex P 



" 3fi ?A^+i4 r ) 



a 2 (ig z )e l «^- J « a!£ exp 



3 ,2 



(69) 



where for z 6 M+ we have 0] 



i 



*K 2 



b{z) := 2^T 



^(A 2 +z) 



/ cosh-y/5 

y'z — tanh -y/z 
P/3 



(70) 
(71) 



For d = 3, the additional degrees of freedom associated 
with excursions in the y-direction lead to the replacement 



of q z a x (k) by q z [a x (k) + al(k)\ which results in an addi- 



tional factor of yl + iq z \ k for each mode k. Thus, for 
general d, we have to replace 02(2) with 



a d [z) := Yl 



1 



1 + zA- 2 



(d-l)/2 



(72) 



As $l[b(iq z )] > for all q z £ [—00,00], the Gaussian 
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integration over q x in Ea. (|69|l can be performed by com- 
pleting the square, such that 



Pd(x,p) 



2tt 



^6irb(iq z 



: exp 



6b(iq z ) 



(73) 



Along similar lines one may also calculate the full joint 
distribution function for a grafted polymer in d = 3, 



P${x,y,z) 



1 



Li in 



^7~ e fi u = \ exp 



x 2 + y 2 



-h{x,y,p) . 



Qb{\q z ) 
(74) 



In addition to the poles of 03(15^) at q z = 1A 2 , on the posi- 
tive imaginary axis of the g z -plane, the integrand also has 
singularities at the zeros 1A 2 , of b(z). Thus we continue 
by evaluating the integrals numerically. 

1. Numerical evaluation of integrals 



leads to a x .t = a/Xk where a is some number. We thus 
find 



and 



resulting in 



As L 2 ± /LLn 



La]T A 



fc=i 



2 _ T " 

~^2 



fc=l 



1 /x* 

2 IT 



1/3, this is equivalent to 
Hence P(x,p) has to vanish for p < x 2 /6 



(77) 



(78) 



(79) 



(80) 



The integrand of Eq. I|73l) has no singularities on the 
real g 2 -axis. Before attempting a numerical integration, 
we discuss the behavior of the different terms appearing 
in Eq. (J73J). For d — 3, we have 



1 



^7 2 



cosh 1 



(75) 



For q z £ K, the real and the imaginary part 1/ cosh-v/i^T 
are respectively even and odd functions rapidly decaying 
in magnitude for q z — > ±00. The real part of l/b(iq z ) 
is strictly positive and increasing with increasing \q z \. 
The imaginary part of l/b(iq z ) behaves asymptotically as 
3 [6~ 1 (iq z )] ~ q z leading to a second strongly oscillating 
contribution to the integrand of Eq. lj75l) besides exp(iqz). 
In the interest of numerical stability of the integration, 
it is advantageous to rewrite the integrand appearing in 
Eq. O to 



1 



3 ig(p-£ 2 /6) _ 



03(19) 



: exp 



2tt ^2^36(1,?) 
for q larger than some fixed go 



- Y (l/b(xq)-xq) 



(76) 



2. Region of vanishing probability 

Eq. H76fl suggests that p = x 2 /6 is a special situation. 
The probability density P(x, z) must vanish for points 
which are at distances greater than L from the graft: 
x 2 +z 2 < L 2 . What does this translate to in the harmonic 
approximation? The largest value £* of r x (L) that can 
be obtained for a given value z* of r z (L) can be found 
from the variation of r x (L) — p(z* — r z (L)) where p is a 
Lagrange multiplier. Using Eq. (|18f) and Eq. (|62|l . this 



3. Results for the general distribution function 

It is now straightforward to evaluate the integrals in 
Eq. H73|) by some standard numerical method. The cor- 
responding results are shown in FigEl as contour plots 
of P(x, p) in d = 3 and d = 2, respectively. These 




FIG. 9: Contour plot of the probability density P(x, p) in 
d = 3 (top) and d = 2 (bottom). 



analytical results compare very well with MC results for 
polymers with a stiffness parameter e < 0.2; see a plot 
with e = 0.1 in FigEH There are deviations between the 
harmonic approximation and MC data for larger values 
of £ HHllL 
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FIG. 10: Color coded density plot of the probability density 
P(x, p) in d = 3 (top) and d = 2 (bottom) obtained from MC 
simulations for e = 0.1. 



The density distribution essentially vanishes outside 
the parabola given by p = x 2 /6, corresponding to the 
classical contour of the polymer in harmonic order. The 
main weight of P(x,p) is concentrated close to this line, 
where the effect is stronger for d = 2. Profiles parallel 
to the p direction are of a shape qualitatively similar to 
P\\(p) (see Fig|2J) at least for small i. Profiles parallel 
to the x-axis are not Gaussian. For small p ;$ 0.1, they 
are peaked at x = but unlike a Gaussian, they vanish 
for x 2 > 6/5. For larger p, they display a double-peaked 
shape. Both features would be completely missed by a 
factorization approximation P(x 1 p) — P±(x)P\\(p). An 
elaborate discussion of the features of the distribution 
function in d = 2 and d = 3 as one increases the stiffness 
parameters or introduces some backbone elasticity will 
be the topic of a forthcoming publication [2(j . 

The shape of the full joint probability distribution 
Ps(x,y, z) is best illustrated by plotting an isosurface, 
e.g. P 3 (x,y,p) — 0.1 as shown in Fig^J Due to ro- 
tational symmetry a density plot for P 3 (x,y,z) may be 
obtained by rotating the contour plot of P%{x, z) (Fig|3J) 
around the z-axis. Again MC and analytical results are 
identical for small e. 




FIG. 11: 3d probability isosurface Ps(x,y,p) — 0.1 
grafted polymer calculated numerically from Eq. 1741 . 



of a 



p are measuring the transverse displacement of the tip x 
and the stored length L — z in units of the characteristic 
transverse and longitudinal length scales, L± and re- 
spectively. As can be inferred from Fig|S]the wall crosses 
the x- and p-axis at 

Lcos-d — ( Lcos-d — ( , . 

Lismp 11 L||coso 

respectively; see FigEl These are the two basic dimen- 




FIG. 12: Geometry of the problem in terms of the reduced co- 
ordinates x and p. The position of the wall is characterized by 
its slope p — r/^/ri± and m, the distance from the origin along 
the S-axis, i.e. the minimal reduced stored length imposed by 
the presence of the constraining wall. To harmonic order the 
finite length of the filament also constrains the reduced stored 
length p to be larger than i 2 /6. 



B. Entropic forces: scaling functions 

We are now in a position to evaluate the general ex- 
pression Eq. (|60() for the restricted partition sum. Before 
going into the details of the calculation it is instructive 
to have a look at the geometry of the problem in terms 
of the dimensionless variables x and p. Recall that x and 



sionless variables characterizing the entropic forces ex- 
erted on the inclined wall. We also introduce the slope 
H = tan a of the constraining wall with respect to the 
x-axis 



p = — = tant9 = j_ tan J . 



3s 



(82) 
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As discussed above the finite length of the polymer gives 
a constraint on the reduced stored length p such that it 
has to be larger than x 2 /6, i.e. above the parabola drawn 
in Fig ll2l Hence, just the points on the constraining 
wall inside the parabola are accessible to the tip of the 
polymer. As one moves the wall further away from the 
grafted end, the number of contact points decreases and 
finally reduces to zero when the wall becomes tangent to 
the parabola. In this limit, where 



(83) 



the force exerted on the wall vanishes. 

We may now write the restricted partition sum in 
terms of the reduced stored length m and the slope of 
the wall /i 



where 

Z{ri\\.ji) = 

dq 
o *q 



1 r m 

-eric _ 

2 V2p 



(84) 



(85) 



(a 3 (iq)e 2 



(p,q) 2 b( lq ) _ e -\( m f 



as shown in Appendix [D] The force is again found by 
taking the derivative of k^T In Z with respect to C- It 
obeys the scaling law 



f(t,#,L,£ p ) = f D {0)f(r, h i i ) , 
with an amplitude 

tt 2 k f c 



and a scaling function 



4L 2 cos d cos § 



(86) 



(87) 



4 z'(n hf i) 
n2 z(v\\,ri 



that can be expressed in terms of the restricted partition 
sum and its derivative 



—$t(e iqri Ua 3 (iq)e- 

7T V 



(89) 



As detailed in Appendix [D] Eq. iJSSJ and Eq. |JSSJ are 

suited best for a numerical evaluation of the entropic 
force. 

In FigEHthc analytical results for the scaling function 
f(r)\\,fj,) of the entropic force are shown as a function of 
<5?7|l = ?/|| — rjp for a series of values for /x. Since we have 
subtracted off the critical value of the reduced stored 
length r)p the forces vanish for 6r)\\ < 0. There is a dra- 
matic difference in the shape of the force-distance curves 
in 2d and 3d. Whereas the force increases monotonically 
with increasing 5m for 3d, it shows a pronounced max- 
imum in 2d, the physical origin of which is the same as 
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FIG. 13: Scaling function of the entropic force / in d = 3 (top) 
and d = 2 (bottom) as a function of 6r)\\ = r)\\ + 1/^ 2 for a series 
of values of fi as indicated in the graphs. Solid lines represent 
analytical results as obtained from a numerical evaluation of 
Eq. l88t . Monte Carlo data for a stiffness parameter e — 0.1 
are given as symbols as indicated in the graphs. For [i = 
one recovers the results for /||(?7||) as discussed in Section UTTl 



for = 0. The maximum in 2d vanishes upon increasing 
fi, which can either be understood as an increase in the 
inclination angle or an increase in the persistence length; 
see Eq. 11551. 

For comparison MC data are given for a particular 
value of the stiffness parameter, e = 0.1. In this stiff 
regime the analytical results compare very well with the 
MC data, except for large values in the stored length 
where the harmonic approximation is expected to become 
invalid. 

For small values of \x the reduced stored length m is no 
longer a good variable. Instead we define a new scaling 
function f(rj±, /x) such that 



/(C,0) 



L± sin & 



(90) 
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where 



where 



f(v±,v) = f{v\\Iv,v) ■ 



(91) 



Like in the previous scaling plot the force should van- 
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FIG. 14: Scaling function f(r]±,ij,) in d = 3 (top) and d — 2 
(bottom) for a series of values for fi (solid lines). For large 
(j,, the scaling function /(r/j_,£i) asymptotically converges to 
f±(r]±) obtained within a factorization approximation. The 
MC data indicated by different symbols in the graphs are 
given for a fixed stiffness parameter e = 0.1. 



ish for 8m < 0, which in terms of r)± reads r]± < — 
Again, there is a marked difference between 2d and 3d 
results; see Fig^] We also observe that the scaling func- 
tion f{rj±,fi>) asymptotically approaches a limiting curve 
for fi — > oo, which for a fixed value of e corresponds 
to i9 — ► 7r/2. It turns out, as we will show now, that 
this limiting behavior can well be explained within a fac- 
torization approximation P(x,z) ~ Pu (z)P±(x). Then, 
Z(C,i?) simplifies to 



sin 



-i?-zcot#). (92) 



2 x (ar) = / dx'Px(a:') 

•/ — OO 



(93) 



is the restricted partition sum for the transverse fluctua- 
tions. The longitudinal distribution function Pu (z) is, for 
small L/lp, strongly peaked at z sa L with a character- 
istic width of L||, and Z± varies on the scale L±. Then, 
for fi ^> 1, the width of the longitudinal distribution 
function is much smaller than the transverse restricted 
partition sum, such that the integration over P^ can 

be approximated by Z(£,$) Zj_ ([£ — Lcosfi] sin -1 #) 
which upon using that the transverse distribution func- 
tion is a simple Gaussian, Ea. (|64(l . results is 



Z(C,0)«ierfc^|=: Z X ( V± ) 



(94) 



This approximation fails when fi w 1, which defines an 
angle 



# c = arctan(L||/Lj_) w ^JsL/£ p 



(95) 



well above which the factorization approximation is valid. 
The entropic force is then 



/(C,tf) 



k B T 
L± sin 



where 



f±(v±) = - 



-nl/2 



Z ±(V±) V 7T erfc(ryj_/v / 2) 



(96) 



(97) 



This result for the scaling function of the entropic force 
is indicated as the thick solid line in Fig^] It becomes 
exact in the limit $ = 7r/2, where starting from Eq. (|(j(jp> 
one can integrate out the longitudinal coordinate to end 
up with 



(98) 



Finally, for large £, one recovers the linear response result 

/(C,7r/2)=3</L 3 . 

If we compare the results of the factorization approxi- 
mation for ■& > i9 c , Eq. I® and Eq. 0[JJ, to Eq.(2) and 
Eq.(5) of Ref . [2l| . one realizes that they are almost iden- 
tical up to the minor difference that Mogilner and Oster 
define their ko to be 4k/L 3 where it actually should be 
3k/ L 3 . The factor 4 in Ref.|2lJ instead of the correct 
value 3 is the result of assuming that the minimal en- 
ergy configuration of a thin rod bent by application of a 
force to its non-grafted end has constant radius of cur- 
vature for small deflections, which is not the case. In 
fact, the boundary condition of the mechanical problem 
forces the curvature to vanish at the non-grafted end. In 
Ref. [IJ the entropic force was calculated by taking into 
account transverse fluctuations of the grafted polymer 
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only and completely disregarding any stored length fluc- 
tuations. Here, the factorization approximation, which 
treats longitudinal and transverse fluctuations as inde- 
pendent, gives the same result for inclination angles 
i? > i? c . The reason behind the validity of the asymptotic 
results, Eq. (|96(l and Eq. Ij97(l . is that the tip distribution 
function is much narrower in the longitudional than the 
transverse direction for $ ^> d c ~ y 'L/£ p . Hence the 
range of validity of the factorization approximation be- 
comes larger as the polymers become stiffer. Of course, 
the analysis in Ref.|2l[ has to fail for small inclination 
angles since it does not account for stored length fluctu- 
ations at all. This is seen most dramatically for d = 0, 
where such an approximation would give no force at all 
in contrast to what we find in Section IlIII 

C. Entropic forces: explicit results 

The analysis in the previous section gives the full scal- 
ing picture for the entropic forces as a function of the 
scaling variables rj\\ and Here we discuss our find- 
ings in terms of the actual distance of the grafted end 
to the wall £ 5 the inclination angle d, and the stiffness 
parameter s = L/£ pi which may be more convenient for 
actual applications. Of course, the disatvantage of such 
a representation is that we now have to give the results 
for particular values of the stiffness parameter. In this 
section we also restrict ourselves to the discussion of fil- 
aments which are allowed to fluctuated in 3d. 

In FigsEl an d El the force / in units of the Euler 
buckling force f c is shown as a function of £ (in units 
of the total filament length L) for a series of values of 
d and vice versa; the stiffness parameter has been taken 
as e = 0.1. Recall that the angle d = corresponds 
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FIG. 15: Analytical and MC simulation results for the en- 
tropic force f/fc as a function of the distance of the grafted 
end from the wall C,jL for a series of inclination angles 
■& — 17, ■ ■ ■ ,89 with steps 9 (in degree). 



to a wall perpendicular to the orientation of the grafted 
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FIG. 16: Analytical and MC simulation results for the en- 
tropic force f/fc as a function of the inclination angle i? 
(in degrees) for a series of distances to the wall C,/L = 
0.99,0.985, ■•■ ,0.95. 



end of the polymer, which has been discussed in detail in 
Section Hill Upon increasing the inclination angle •& the 
entropic force decreases for all given values of (. This is 
to be expected since the wall then cuts off less from the 
probability cloud of the polymer tip. For the same reason 
the forces also decrease with increasing £ for a given value 
of t9. The analytical results (solid lines) agree well with 
the MC data for not too small values of £. The deviations 
grow larger upon decreasing the distance between the 
wall and the grafted end. Then non-linear effects not 
taken into account by our weakly bending approximation 
set in. 

In the limit as the inclination angle approaches ir/2 it 
is certainly no longer justified to calculate the entropic 
force by assuming that only the polymer tip is not al- 
lowed to penetrate the membrane. Then, one has to take 
into account the fact that also the body of the polymer 
is constrained by the presence of the wall. Since this re- 
duces the number of allowed polymer configurations even 
further this effect is expected to lead to an enhancement 
of the entropic force. Indeed this is the case, as one may 
infer from Fig ll7l where we show a comparison with MC 
simulation accounting for these constraints. One also 
notes that the enhancement of the entropic forces be- 
comes largest as t9 — > ir/2 and the distance between the 
wall and the grafted end becomes small; a full account of 
this effect will appear in Ref. 22] . 

Finally, we would like to compare our full results with 
the factorization approximation discussed in the previous 
section, Ea. (|97(l . which when corrected for some minor 
factor is identical to the results given in Ref. [21] . The 
comparison is given in Fig llHI for a stiffness parameter 
e = 0.1. In the limit of large inclination angles well 
above -d c ss 30°, there is excellent agreement between the 
factorization approximation and the full results for not 
too small values of £. As one approaches d c the range of 
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FIG. 17: Comparison of the analytical results for the entropic 
force as a function of £/L (solid lines) for a series of values for 
$ = 17, • ■ ■ ,89 with steps 9 (in degree) indicated in the graph 
with MC simulations (symbols as indicated in the graph), 
which take into account the constraints of the wall on the 
body of the polymer; e = 0.1. The deviations are most pro- 
nounced for small values of ( and inclination angles $ close 
to tt/2. 



□ 
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FIG. 18: Comparison of the full result for the entropic force 
as a function of £ (full lines) with the results obtained from 
the factorization approximation (dashed lines) for a stiffness 
parameter e = 0.1 and a series of inclination angles $ = 
30° , 40° , • ■ • , 80° . The range of validity of the factorization 
approximation broadens as one increases the inclination an- 
gle #; it is invalid for ■& smaller than # c ~ 30°. 



validity of the factorization approximation shrinks and 
finally it becomes invalid for d < d c . 

To illustrate the applicability of the factorization ap- 
proximation let us take some examples. For the cy- 
toskelctal filament F-actin with a contour length 100 nm 
and persistence length £ p = 15/^m, the stiffness param- 
eter becomes e — 0.006 which gives i? c w 7.6°. Upon 
increasing the stiffness parameter to £ = 0.1, which 



amounts to changing the contour length to a value of 
L = 1.6/im, the critical angle # c increases to 28.7°. 



V. SUMMARY AND CONCLUSIONS 

In summary, we have presented analytical calculations 
and extensive Monte Carlo simulations for the entropic 
force / exerted by a grafted polymer on a rigid obstacle 
(wall) . The scale for the magnitude of the entropic force 
is given by the Euler buckling force f c cx k^Ti v / L 2 . The 
stiffness parameter e = L/£ p discerns the two universal 
regimes of a Gaussian chain (e » 1) and a semiflexible 
chain (e < 1). In this manuscript we have mainly fo- 
cused on the stiff limit, where analytical calculations us- 
ing a weakly bending rod approximation are possible. In 
comparing our results with Monte Carlo simulations we 
have found that the range of applicability of the results 
obtained in the stiff limit extend to stiffness parameters 
as large as e = 0.1. Qualitatively the asymptotic results 
remain valid even up to e = 1. 

For the simplest possible geometry, where the polymer 
is perpendicular to the wall, located at a distance £ from 
the grafted end, our analytical calculations show that the 
entropic force obeys a scaling law in the stiff limit 



/iiCCMpWc/u^) 



(99) 



with the scaling variable fj = (L — 0/L\\ measuring the 
minimal compression of the filament in units of the longi- 
tudinal width of the tip distribution function Lu = L 2 /£ p , 
and f c the Euler buckling force of a classical beam. For 
small values of the scaling variable we have derived a 
simple analytical expression, Eq. (14-51) . 



fHv) 



exp[— 1/477] 



7,5/2 ^3/2 [1 _ 2erfc( 1/VrJ)] 



(100) 



and a corresponding formula in 2d, Eq. (|58|) . which de- 
scribe the full scaling function to a high numerical accu- 
racy for fj < 0.2. For fj > 0.2 there are equally simple 
expressions, as for example Eq. (|46fl for 3d. We expect 
these formulas to be useful for molecular models of cell 
motility. The shape of the scaling function shows dra- 
matic differences between 2d and 3d, which are of geo- 
metric origin. In 3d the entropic forces always stay below 
the Euler buckling force. In contrast, in 2d it is larger 
than the mechanical limit for most of the parameter space 
and exhibits a pronounced maximum at small values of 
the scaling variable fj before it steeply drops to zero as 
C - L. 

Extensive Monte Carlo simulations confirm these ana- 
lytical results and show that their range of applicability is 
e < 0.1. For larger values of the stiffness parameter there 
are clear deviations from the stiff scaling limit, which be- 
come qualitative for e > 1. Features of the stiff limit, 
such as the maximum in the entropic force, are visible 
even for e as large as 4. 
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For a polymer inclined at an angle $ with respect to 
the wall also the transverse width L± = <JL 3 /3£ P of the 
tip distribution function plays a significant role; note that 
the ratio L\\/L± — \/3e. The entropic force can now be 
written in the scaling form 

f(C,#;L,£ p )=f c (#)f(r) h r) ± ), (101) 

where r/± = (L cos d — Q/(L± sinz?), = (Lcost) — 
£)/(L|| cost?) and / c (i?) = f c / cosi3. It turned out that 
a proper choice of scaling variables are fi = m/rjx = 
(L^/Ln) tarn? and 7711 or rj± depending on whether the 
inclination angle is smaller or larger than a characteristic 
angle tani9 c = Lu/L±, i.e. fi c — 1. Upon increasing the 
inclination parameter [i the shape of the scaling function 
changes from a step-function-likc form at ji = to a 
purely convex shape as fi — ► 00. The limit fi — ► 00 cither 
corresponds to 1) — > tt/2 or for a fixed 1? 7^ to the stiff 
limit £ — > 0. For 2d, in addition, the maximum vanishes 
at w 0.6. 

In the limit of inclination angles which are much larger 
than the characteristic angle $ c , we have found that an 
approximation, Eq. (|96|l and Eq. (|97|l . based on factoriz- 
ing the joint probability distribution of the polymer tip 
gives an excellent asymptotic representation of the full 
analytical results: 

^-LlS^n^m- (1 ° 2 ' 

It is simpler than the full scaling form since it only de- 
pends on one scaling variable. Up to minor factors this 
asymptotic formula for the entropic force is mathemati- 
cally identical to the results found in Ref . [2l| , which was 
derived upon assuming that the tip of the polymer fluctu- 
ates perpendicular to its contour only. Since tan i? c oc *J~e 
the range of applicability of this results grows with in- 
creasing stiffness parameter. For example, d c equals ap- 
proximately 30° and 10° for stiffness parameter e equal to 
0.1 and 0.006, respectively. For d < $ c the factorization 
approximation fails completely, since it gives an incorrect 
description of the longitudinal stored length fluctuations. 



Then, a full analysis in terms of a two parameter scaling 
function is necessary. 
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APPENDIX A: INVERSE LAPLACE 
TRANSFORM OF THE MOMENT GENERATING 
FUNCTION 

In this appendix we collect our calculations of the in- 
verse Laplace transform of the moment generating func- 
tions. This will give as two sets of series representations, 
which show good convergence properties either close to 
full stretching or for strong compression of the filament. 



1. Series representation of the 3d tip distribution 
function for large stored length 

Starting from the moment generating function V\\ (/) 
one can calculate the distribution function Py (z) by an 
inverse Laplace, i.e. an integral along the imaginary axis, 

/~\~ 1 CO 1 P 
^ ^: e/(L ~ z)7 W)- ( A1 ) 

Since the moment generating function 

has poles at fk = — \\l v jL 2 with k = 1, 2, 3, ... only along 
the negative real axis, standard residuum calculus gives 



P||(/)=£exp 



fe=i 



<L-')\lf, 



n 1 



l^k 



(2k -If 
(21 -If 



f A 2 



using nr =1 (1 



cos (-fx) [T3 |. the product term can be written as 



n 



l^k 



(2k- l) 2 

(21 — l) 2 



»mii-^)n(i (2t '" 1)2 



k'^k 



(2k- 



, / (2k 1 - I)' 
hm 1 — ; -7- I cos 

k'^k \ (2k- l) 2 

2 (-l)* +1 2 =2( _ 1)fc+1 1 
7T 2k -I y ' Xk 



(21 I) 2 

1 (^-1 



(A3) 



(A4) 
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Hence we find 



P|| (z) = 2L- 1 £(-1)*+% cxp [-X 2 k (L - 2)/L,|] 



(A5) 



fc=i 



with the characteristic longitudinal length scale L\\ = L 2 /£ p . 



2. Series representation for the tip distribution function close to full stretching: general d 



We begin the analysis with the two-dimensional case, where 



k=l 



fL 



-1/2 



1 



A fc J ]l cosh fL\\ 



(A6) 



For the derivation of our first series representation we start from the product formula for the moment generating 
function. In this representation one has branch cuts on the negative real axis at / — fL\\ = —X 2 for k E N. We now 
deform the contour in the complex plane such that we enclose the negative real axis. Then 

' 00 df 



P\\(P) 



2m 



df 



2?ri 



^-e"P|,(/-ie) + \ -L.^"V\\{f + ie) 
2tti 



-fp 



(A7) 



where e — + 0. To proceed we need to evaluate the product formula on the negative real axis. We find for x G 
[2fc + l,2fc + 3]§ 



i=i 



x 2 =F le 



(t0* 



i 



^/| cos: 



Upon substituting y 2 = f this finally results in the series expansion 

2 00 ^2 
7T 



^ J\ 2n + i VI COST/I 



71=0 



ye 



(A8) 



(A9) 



For large values of p, corresponding to a significant compression of the polymer, the integral is dominated by the 
contribution from the interval [tt/2, 3tt/2], such that the leading factor will be proportional to exp(— 7r 2 j 5/4). In order 
to evaluate P|| (p) further, we may average ye~ v p over the interval and approximate the integral as 



ye~ v p 



^271 + 2 

dy- 

a 2 „+i V\ cos y\ 



+ » exp 



m=4 



*2n+? P 



dy 



\/\cosy\ 



(A10) 



such that we finally get 



-A 



2>H 



n=0 



where 



07T 



V|cosy| 
I 



0.67. 



(All) 



(A12) 



Next we drive a series representation suitable for small With cosh(x) = ^(e^+e x ) and the generalized binomial 



values of p. We use that for / e K + one has 01 

W) 1 



theorem, this can be expanded to give 



cosh yf]L\\ 



(A13) 



1=0 



2 I e -(2'+l/2)V/ i ll 



(A14) 
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which is a holomorphic function on C\M_. Hence by 
the theorem of identity from complex calculus this for- 
mula remains valid V/ S C\R_. Substituting y = y/fL\\ 
transforms Eq. IjAlfl to 



P \\(P) 



dy 



e y ?y-P\\(y 2 )- 



(A15) 



— lOC + € 



Inserting the series representation Eq. (| A14|) and using 
the integral representation 



Di(z) = V2ttc 



ioc+e 



ioo+e 



ds 
2tti 



s exp 



for the parabolic cylinder function [16j as well as 

(S) =(-!)' (2/ " 1)!! 
where n\\ — n(n — 2)(n — 4) . . . yields 



(A16) 



(A17) 



x exp 



2 l l\ 





Di 




[ 2p _ 




[ V2~P \ 



(A18) 



With Di(x) = xe-* 2 ^ Eq. llAT8)l becomes Eq. (|5U|) . 

Finally, all the calculations are easily generalized to 
general spatial dimensions d. One finds the series repre- 
sentation 



OO / 1/7 \ 



(=0 



J + i(d-l) 



p3/2 



exp 



(I 



(A19) 



which is the fast converging for small p. 

APPENDIX B: SADDLE POINT 
APPROXIMATION 

Starting from Eq. IjAlf) and introducing / = fL^ gives 

/ + 10C ir 
_ i^cosh-^ 

(Bl) 



ioo 27rl eV? + e~V? 



We are interested to the asymptotic result of the integral 
close to full stretching p — > 0. Upon substituting / = 
£//5 2 one finds 



W 4n (B2) 



where /(£) = £ — v 7 ?- Since the function /(£) has a 
global maximum at £o = 0.25 the main contribution 
to the integral in the limit l/p — ► oo comes from the 
integration along the curve of steepest descent which 
passes through £q- We need to find this curve such 
that $[/(£)] = constant = $*[/(&)] = 0. We write 
v 7 ^ = y^a(l+ is) in terms of the curve parameter s. Then 
the condition 3[/(Co)] = gives a — 1/4, and the curve 
of steepest descent is given in terms of Sft[£] = |(1 — s 2 ) 
and 9f[£] = 2as, which is a parabola parameterized by 
s. The saddle point approximation amounts to a contour 
integral along this parabola, where /(£) = — (1 + s 2 )/4, 
such that 



1 



L \\P 2 



ds e -(i+s 2 )/±P 

^+"> 1+e - w „)/ 8 - < B3 > 



To the leading order in p we get 



W = 



exp[-l/4p] 



P 2L \\ 



ds 

27 6XP 



■ exp 



1 



(B4) 



In the two dimensional case (2d), using the same strat- 
egy and substituting / = ^/p 4 ^ 3 gives 



exp 



1 

16/5 



(B5) 



APPENDIX C: JACOBI TRANSFORMATION OF 
THE RESTRICTED PARTITION SUM Z\\ (0 

To unclutter the formulas in this section, we use the 
generic argument x with x = rju . Z\\ (x) can be written 
as 

/■OO °° 



Z\\{x) =2 dy V (-l) k+1 6(y - \ k )- e -**v 2 * 

J " k =-oo y 

' dyS(y)-e- e -y 2x (CI) 
V 

where we defined 

oo 

6(3,):= J2 (-l) k+1 S(Xk -y). (C2) 



k=—oo 



Since S(y) is odd in y and has periodicity 2tt, we can 
expand it into a Fourier-sine-series: 



Hy) =^2dism(ly) 



(C3) 
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where 



This results in 



di = - dyS(y)sm(ly) 
n Jo 
2 

= -sin(77r/2) 

2 J if Ms even 

7r 1 (— if l i s odd 



(C4) 



APPENDIX D: GRAFT-ANGLE-DEPENDENT 
FORCE 

We evaluate the general expression Eq. i|60J) using the 
representation 



0(x) = lim / 

6^0+ J 



dq e lqx 



2m q — ie 

of the step function Q(x). With Eq. (O we find 



(Dl) 



*(i/) = -£(-i) I+1 siii[(2i-i)i,]. 



(C5) 



Inserting this into Eq. (|C1(I we find for Z\\ (a 



A °° />0O 

J ll( a; ) = -E(- 1 )' +1 / ^y'e-^sin^-l)^ 
71 i=i ^° 



(C6) 



The integral evaluates to 0] (with /z = 0, /3 = a;, 7 
2Z— 11 



dj/l/-^"" *sin[(2Z-l)y] 

(2Z - i) e -(2/-i) 2 /4^ / 3 (21-1)2 
2~T X ^ lFl [ 1] 2 5 ~^x~ 



(C7) 



As the confluent hypergeometric function iiq(a;7;z) = 
<j>(a , 7; 2) has the property &(a; 7; z) — e z $(7~ a, 7; — z) 
|16| we find with [13 



rfv^. 



2^1 



Our result for Zn (x) is thus 



i,|( a; ) = 2£(-l) i+1 erf 



21-1 



(C8) 



(C9) 



This still has problems for x — ► where erf[(2Z 
l)/2-y/x] — > 1. We can, however rewrite it to 



21 - 1 



i,|(x) = 2^(-l)' +1 + 2 £(-l) I erfc— ^ 

i=i ;=i "^^ 



(CIO) 



All convergence problems are now isolated in the first 
sum. As we know that -Zn(O) = 1 (compare Eq. wc 
assign 2^^ 1 (— 1) /+1 = 1 to finally find 



Z\\(x) = l + 2^(-l)'erfc 



1=1 



21-1 

2^X~ 



(Cll) 



dq ex P 



C/cosi?-L 
^Lv 



27ri q — ie 

dxdr)e iqf >e- iq{L ^ /L ^ tanm P{x,fi) 



dq 



exp 



C/cosi9-L 




2tti 



g — ie 



exp 



(qi ± L | [ 1 tan^) 2 36(-i(7)/2 



(D2) 



where 



Z(r, h ti = / ^^a d ( ig )e ^ . (D3) 

J 27ri g + is 

Using the Dirac formula 



_L_ = V l _ 17r( J(g) 

q + ie g 



(D4) 



03(0) = 1, 36(0) = 1 and the symmetry properties of 
03(19) and b(iq), we find 

Z(„||,m) = \- 2 j™ ^ (e^a 3 (iq)e-^ 2q23b ^)(D5) 

The notation "P denoting the principal value has been 
dropped as the integrand is regular at q = 0. For large /z 
and/or £, -Z(?7|| , /i) vanishes. This means that the integral 
in Eq. (|85|l must approach 1/2. Subtracting the result 
of the numerically evaluating the non- vanishing integral 
from 1/2 strongly amplifies the unavoidable round-off er- 
ror. We therefore rewrite Eq. I|D5I) to 

Z(ni,u) = -erfc ^!! 



-2 / -i-9f 
/o 27r 9 



I a 3 (ig)e"^ 2 ' z23fc( " 3) 



e 2 



where we used the identity 
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2- V 



dq e 19?) ii 

— — < 

27T1 Q 



ierfci 
2 y/2fi 



(D6) 



(D7) 



21 



As Qq 2 b(iq) ~ — q for large \q\, it is again advantageous 
to split the integrals at some qo and, for q > q , to rewrite 
the imaginary part appearing in the integrand of Eq. l|D6l) 
to 



Hence the integrals vanish if 6r)\\ := r/|| + 3/i 2 /2 < which 
we already understood in the simple geometric picture of 
the problem. 



Both integrals now vanish in the limit of large rj^ and 
(D8) have well-behaved integrands on [0, oo]. The precision 
with which /(7711,/x) can be calculated is, however, still 
limited by the relative error in evaluating the integrals. 
This relative error grows quickly with increasing m limit- 

(D9) ing the range of rju over which /(1711, /u) can be calculated 
reliably (note that the first term of Eq. (|D6|I vanishes 
In both cases the integrand is holomorphic for 3g < 0. with increasing jjii as well). 



and the real part appearing in Eq. 189(1 to 

K (eMvn +3m 2 /2) as ( iq ) c -^ 2 (9 2 b(>?)+>?) 
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